CERN - European Organization for Nuclear Research 



LCD-Note-2011-003 

Determination of Chargino and Neutralino Masses in 
high-mass SUSY scenarios at CLIC 

N. Alster* f and M. Battaglia* f 

University of California at Santa Cruz, Santa Cruz, CA, USA ' CERN, Geneva, Switzerland, 

March 11,2011 
Abstract 

This note reports the results of a study of the accuracy in the determination of 
chargino and neutralino masses in two high-mass supersymmetric scenarios through 
kinematic endpoints and threshold scans at a multi-TeV e + e— collider. The effects 
of initial state radiation, beamstrahlung and parton energy resolution are studied in 
fully hadronic final states of inclusive SUSY samples. Results obtained at gener- 
ator level are compared to those from fully simulated and reconstructed events for 
selected channels. 



1 Introduction 



The study of the gaugino sector of Supersymmetry is a complex and important endeavour, which 
appears well suited to a linear collider of sufficient energy and luminosity. The main observables 
of interest are the masses of the x° an d X ± states and their production cross sections, including 
those with polarised beams. e + e~ collisions offer two independent techniques for determining 
the mass of supersymmetric particles. These are the analysis of the energy spectrum of the SM 
particle produced in association with a lighter supersymmetric state in the two-body decays and 
the study of the pair production cross section near threshold. These techniques have already been 
extensively studied for lower centre-of-mass energies, y/s, between 0.35 to 0.5 TeV HJHl. m this 
note, we analyse the gaugino pair production and derive the statistical accuracy on their masses 
using both techniques and including the effects of initial state radiation (ISR), beamstrahlung 
(BS) and parton energy resolution for multi-TeV e + e~ collisions. We follow the evolution of 
these accuracies for fully hadronic final states from pure signal samples to realistic inclusive 
SUSY samples and validate the results obtained at generator level with analyses performed on 
fully simulated and reconstructed events. The study provides us with requirements on parton en- 
ergy resolution which are complementary to those obtained from other processes, such as heavy 
SUSY Higgs decays, since the kinematics of decays of gaugino pairs with large missing energy 
into pairs of escaping neutralinos does not benefit from the kinematic fits, which are instead 
applicable to processes where the full beam energy is deposited in the detector. The estimated 
mass accuracies can be compared in a next step to those required for the reconstruction of the 
GUT scale SUSY parameters and the determination of the lightest neutralino contribution 
to the dark matter relic density in the universe [6 |. This comparison will provide us with well- 
motivated quantitative requirements on parton energy resolution in SUSY events. 

2 SUSY Models 

This study considers two scenarios in the constrained MSSM (cMSSM) model , which offer 
different experimental challenges. Their parameters are given in Table [T] The first (model 
I), adopted as a benchmark point for the CLIC CDR studies Q, has the lightest neutralino at 
340 GeV and the chargino and heavier neutralinos with masses in the range 640 to 917 GeV 
(see Table [2] and the left panel of Figured}. At \/s = 3 TeV all the gauginos are observables. 
The relatively low masses and the 3 TeV centre-of-mass energy make cross sections sizable but 
the beamstrahlung effects more significant (see Table[T}- In the second (model 110) the lightest 
neutralino has a mass of 554 GeV, while the other neutralinos and the charginos have masses in 
the range from 1064 to 1414 GeV (see Tableland the right panel of Figured]) IM- At 3 TeV, most 
gauginos are close to threshold for pair production and cross sections are small. This minimises 
the beamstrahlung effects, since the production cross section drops significantly when the beams 
lose energy due to radiation. The cross sections are given in Table [3] and Figure |2] 

The main features of these two models of interest in this study are the large sparticle masses 
and the dominance of decays through W , h and, to a lesser extent, Z° bosons. These fea- 
tures are common to most of the large-tan /3 cMSSM models (9) with neutralino dark matter 

'This is point K' of ref (8). 
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Table 1: Parameters of the two cMSSM models adopted in this study 



Parameter 


Model I 


Model II 


m (GeV) 


966 


1001 


m 1/2 (GeV) 


800 


1300 


tan/3 


51 


46 


A 


0. 


0. 


sgn(At) 


+ 




m top (GeV) 


173.3 


175 



Table 2: Gaugino mass spectrum in the two cMSSM models adopted in this study 



Particle 


Mass 


Width 


Mass 


Width 




(GeV) 


(GeV) 


(GeV) 


(GeV) 


X\ 


340.3 




554.3 




X°2 


643.2 


0.02 


1064.2 


0.04 


i 


905.5 


4.55 


1407.2 


6.75 


X4 


916.7 


4.64 


1413.8 


6.85 


xt 


643.2 


0.02 


1064.3 


0.04 


X? 


916.7 


4.63 


1413.7 


8.08 



compatible with the Cosmic Microwave Background (CMB) data. We verify this by perform- 
ing scans of the cMSSM parameter space around both benchmark points to examine the mass 
spectrum and decay branching fractions of gauginos. In these scans we vary mo and m\n within 
±300 GeV from the benchmark parameters, tan/3 within +j 5 , Aq within ±250 and \l to have 
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(a) Model I (b) Model II 

Figure 1: Supersymmetric particle spectra for Model I and II. 
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Table 3: Cross sections for gaugino pair production in the two cMSSM models adopted in this 
study 



Process 


no Rad 


ISR 


ISR+BS 


Pol +0.8/0.0 


Pol +0.8/-0.6 




(fb) 


(fb) 


(fb) 


(fb) 


(fb) 


Model I 












Inclusive SUSY 


103.3 


97.4 


79.3 






e + e~ -> XtXi 


11.5 


11.6 


11.9 


21.4 


34.0 


e + e- -> X2X2 


4.2 


4.1 


3.6 


6.5 


10.4 


e + e- -> Z 2 +^ 


14.5 


14.4 


13.8 


21.1 


33.3 




6.1 


6.1 


6.0 


6.8 


10.3 


Model II 












Inclusive SUSY 


45.4 


39.0 


25.2 






e + e- -> jfi^f 


4.9 


4.3 


2.9 


5.2 


8.3 


->■ Z2Z2 


1.9 


1.7 


1.1 


2.0 


3.1 




6.2 


4.6 


2.1 


3.2 


5.1 




2.6 
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(a) Model I (b) Model II 

Figure 2: Cross sections for gaugino pair production vs. y/s in (left) model I and (right) model II 



either sign. The sparticle spectrum corresponding to each set of parameters is computed us- 
ing SuSPECT 2.2 (TOj and the decay branching ratios with SDECAY 1.3 ifTTl . We compute 
the neutralino relic density, Q. x , using Micromegas 2.2 |[T2l . In total, we generate over 150k 
cMSSM points using a flat sampling of the parameter phase space. Of these, we retain those 
points consistent with the current limits on the lightest Higgs bosons and sparticle masses and 
yielding 0.093< Q.%/h 2 <0. 129, in agreement with the WMAP seven-years data [13]. Figure[3] 
shows the % 1 mass and Figure [4] the %i ~^ W ± X\ branching fraction for the accepted points. 
These results show that decays into bosons are dominant in this region of the cMSSM param- 
eter space. These results can be extended to the general MSSM where a significant fraction of 
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Figure 3: Distribution of %f~ mass values for points compatible with WMAP data from the 
cMSSM scan around (left) model I and (right) model II. 
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Figure 4: Distribution of branching fraction values for (left) %\ W ± Xi and (right) x\ ~~ ^ 

Z°Xi for points compatible with WMAP data from the cMSSM scan around (upper 
row) model I and (lower row) model II, showing the predominance of decays into 
bosons. 



the decays of %f and %® and the majority of those of the heavier chargino and neutralinos are 
two-body processes with emission of a boson. In these scenarios, W , Z° and h production is 
a distinctive signature of gaugino decays. Figure |5]shows the boson mass spectrum at generator 
level in inclusive SUSY events for the two models considered in this study, showing the ± W, 
Z° and h° boson contribution. The W/Z/h discrimination is essential for identifying the decay 
processes, which sets constraints on parton energy resolution through the di-jet invariant mass 
resolution. Precise parton energy reconstruction is similarly required to preserve the accuracy in 
the gaugino mass measurements through the determination of the endpoints of the boson energy 
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Figure 5: Boson spectrum in inclusive SUSY events with > 4 jets + missing energy for model I 
(continuous line) and II (dashed line). The h mass is 118.5 and 122.9 GeV, respec- 
tively. 

spectra. 

2.1 Event Simulation 

For this study, events are generated with Pythia 6.215 |[T4l . For model I the SUSY spectrum 
and the sparticle decay branching fractions are forced to those obtained with SuSpect 2.2 and 
SDECAY 1.3, respectively. For model II the spectrum is generated with ISASUGRA 7.69 US- 
Polarised cross sections are obtained using SUSYGEN 3.0 lfl6l . The production cross sections 
for unpolarised and polarised beams are summarised in Table [3] Samples of inclusive SUSY 
events are processed through full detector simulation using the GEANT-4-based MOKKA ifrTl 
program and reconstructed with M ARLlN-based lfT8l processors for a version of the ILD detector 
concept O, modified for physics at CLIC lfl9l . These events are used for a validation of the 
results in the 4-jet, WW and hh final states, discussed in section 3.3.5. 

3 Mass Determination by Fits to Boson Energy Spectra 

In the two body decay process A — > BC into a boson B and a lighter gaugino, C, which are a 
signature of these high-mass benchmark points, the masses of the parent and daughter gauginos 
can be extracted from the position of the kinematic edges of the boson energy spectrum. The 
technique was first proposed for squarks 11201 and later extended to other sparticles in two-body 
decays J2T|. In the case of gaugino decays into bosons, their mass, mj, cannot be neglected, as 
in the case of squark and slepton decays and the relation between the energy endpoint and the 
masses of the particle involved in the decay process are given by: 

E bh , bl = y(E* b ±I3E* b ) (1) 
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where 



m5 + m\ — m 2 r 
A B c (2) 



2ma 



2ma 



(3) 



1^ (4) 

These formulae can be extended in a straightforward way to the case in which the gaugino 
A is not directly produced in the e+e~ collisions but originates from the decay of an heavier 
particle, A', by replacing s with E\, where Ea is its energy. In the case of cascading decays 
A' — > AB' — > BC, Ea is obtained as — E B i H < Ea < y/s — E B i L . 

3.1 Channels 

We study the following processes for model I: 

• e+e- -> XiXi -> W+XyW-Xi; W -> 

• e+e- -> * 2 °* 2 ° -> /i ^^ ^?; A -> bb, 

• e+e" -> ^+^+ -> W+^W-^o -> W+fcO^-^p. A -> W -> w', 

• e+e- -> ^+^+ -> h°xfW-%° -> /i ^^-^ ; A -)• W -> <?<f , 

• e+e" -> ^ 2 +^ 2 + -> Z^+W-^f -> ZOW+^W-^?; Z -> W -> ^, 

• e+e- -> -> W+/f -> W+W-XxW-W+x*; W -> 
For model II we study: 

• e+e- -> -»• W+^OW-^;?; W -> ^ 



e ' e 



z o z o ^ /^p^p. A 



They explore increasing event complexities from 4-jet with two bosons to eight jets with four 
bosons. Mass values and their statistical uncertainties are extracted by a binned x 2 fit using the 
Minuit package l22l . The boson energy spectrum from simulation is compared to reference 
spectra generated according to Eq.(fl3 above, changing the sparticle masses. An integrated lu- 
minosity of 2 ab" 1 is assumed. The y/s value in Eq.© is either kept to the nominal value of 
3 TeV or smeared to account for ISR and beamstrahlung effects, as discussed below. Since the 
fit procedures requires to generate a large number of reference spectra, which are statistically 
independent, in the fitter iterations, each of these is filled with 2.5 x 10 6 random entries dis- 
tributed according to Eq.([T}, to minimise effects from their statistical fluctuations in the Minuit 
calculations, in particular in the determination of the derivatives. 



throughout the paper the charge conjugate of the given state is also implied 
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3.2 Effect of Beam Spectra and Energy Resolution 

The effects of initial state radiation (ISR), beamstrahlung (BS) and finite resolution in parton 
energy reconstruction are taken into account. Beamstrahlung effects on the luminosity spectrum 
are included using results of the CLIC beam simulation for the 2008 accelerator parameters Il23l . 
Initial state radiation is included in the event generation in Pythia. The beamstrahlung spec- 
trum obtained is then used for smearing the y/s value in Eq. 0] We model the ISR spectrum 




Figure 6: Centre-of-mass energy distribution including (a) ISR and (b) ISR and beamstrahlung. 

The points represent the simulation and the lines the phenomenological functions used 
for describing their shape. 




Figure 7: W ± energy spectrum in the % x — > W ± Xi process without radiation effects (dashed), 
with ISR only (dotted) and with both ISR and beamstrahlung (continuous) for (a) 
Model I and (b) Model II. 

by an approximate solution to the Gribov-Lipatov equation, proposed in Il24l . In the formula 
we leave free the rj parameter and the fraction of events off the full energy peak. We determine 
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them by a fit to the ISR spectrum obtained for Pythia signal events (see Figure[6]). We study the 
accuracy on the mass measurements in some of the channels for no radiation effects, ISR only 
and ISR + BS. Figure [7] shows these effects on the energy spectrum for the %t ~^ W ± Xi 
process. In this study we consider only fully hadronic decays of bosons, since we need to recon- 
struct their energy. The effect of the finite resolution in the determination of the boson energy 
is included by applying a Gaussian smearing to the energy of the partons produced in the boson 
decay. The smearing of the energy of the two partons is independent and we test the effect of 
various energy resolution values in the range < 8E/E < 0.125. The energy smearing does not 
account for energy missing in neutrinos. This is particularly important in the reconstruction of 
h° — > bb decays, where either a b or a c hadron decays semi-leptonically. 

3.3 Analysis and Results 

We perform the study by selecting for each of the production and decay processes given above 
the final state parton topologies and boson contents which have the most favourable signal con- 
tribution. We consider only fully hadronic final states and estimate the signal purity and the 
dominant SUSY background processes at generator level assuming perfect jet clustering and 
gauge boson identification. The fraction of SUSY final states yielding a 2-q topology is 4.7% 
(2.3%), a 4-q is 13.3% (1.8%), a 6-q is 3.3% (0.5%) and an %-q is 4.0% (0.8%) for model I 
(II), respectively. In model II, the branching fractions of decays into a^-y t^Xi are impor- 
tant, which explains the reduced rate of multi-quark final states. However, these modes are not 
considered here, since our study focuses on decays into bosons. 

We assume to operate the collider at 3 TeV for an integrated luminosity of 2 ab~ 1 with unpo- 
larised beams. Mass fits are performed on samples of bosons in the selected topology populated 
with either signal only or inclusive SUSY events and we study the evolution of the statisti- 
cal accuracy on the masses with the smearing of the energy spectrum due to beam radiation 
and detector energy resolution effects. We estimate the change in signal purity with the parton 
energy resolution and the identification of the topology using only reconstructed quantities in 
section 3.4. 

3.3.1 %f^W ± x Q i 

The e + e~ — > X\X\ ~^ W + %fW ~X\ ; W — > qc( process leads to a distinctive final state with 
four hadronic jets and missing energy. The SUSY background in model I is almost entirely 
due to charginos produced through e + e~ —?■ e^e^, ef — > %\ v e> which account for 12% of the 
inclusive WW + missing energy SUSY sample. Since the chargino energy in this process is 
lower than the beam energy the W energy is correspondingly shifted to lower values. In the case 
of model II, the signal e + e~ — > X\%\ accounts for 87% of the W + VK~ + missing energy final 
state with background contributions from XzXz an d %z%\ + cc - The main standard model 
irreducible background is due to e + e~ — > W + W~v e v e , which has a cross section of 124 fb. 
However, the W production in this process is forward peaked while it is central in SUSY s- 
channel signal events. Requiring each W boson to be produced within | cos 6\ < 0.85, where 6 
is its polar angle reduces the WWW cross section to 28.7 fb. The energy distribution of W ± 
bosons in WWw events within the angular acceptance, generated with Comphep 4.5.1 IT251 . is 
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E m(55lng (GeV) E m(5slns (GeV) 

Figure 8: Kinematics of e + e~ — > W + W~v e v e and events with |cos$w| < 0.85 obtained with 
Comphep 4.5.1: (upper left) W energy spectrum and (upper right) missing energy 
spectrum. Missing energy distribution for signal e + e~ — > Xi%\ events for (lower left) 
model I and (lower right) model II. 



shown in Figure [8] This background affects only the lower endpoint of the W energy spectrum. 
Further, the bulk of the signal SUSY events can be separated from this background based on 
the event missing energy. The results of the 1-par fit of the %\ mass performed for spectra 
with no radiation, only ISR and also beamstrahlung effects are given in Table |4] Here, we keep 




E w (GeV) E w (GeV) E m (GeV) 



(a) no Rad. (b) ISR Only (c) ISR + BS 

Figure 9: W ± energy spectrum in the four-jet topology, 2 W + missing energy final state for 
X\X\ signal events in model I at 3 TeV with (a) no radiation, (b) only ISR and (c) 
ISR+BS. The points with error bars are the simulation and the line the fitted spectrum. 

the lightest neutralino mass, M o, fixed to its model value, since it should be independently 
determined in the e + e~ — > Mr ~~ ^ M + M X\X\ process and the corresponding e + e~ — > e^e^ 
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Figure 10: W ± energy spectrum in the four jet topology, 2 W + missing energy final state for 
inclusive SUSY events in model I at 3 TeV with only ISR (a) and ISR+BS (b). The 
points with error bars are the simulation, the continuous line the fitted spectrum and 
the dashed line the fitted contribution from X\X\ signal events. 



Table 4: Statistical accuracy on X\ mass from a fit to the E w ± endpoints for X\X\ signal events 
and 2 ab -1 of integrated luminosity at 3 TeV, under different assumptions 



Particle 


Mass 
(GeV) 


NoRad ISR ISR+BS ISR+BS ISR+BS 
8E/E=0 =0.025 =0.05 


Model I 

xt 


643.2 


± 0.91 ± 1.39 ± 2.09 ± 2.89 ± 3.60 


Model II 


1062.2 


±6.10 ±8.25 ±10.11 ±11.0 ±12.0 



to an accuracy of ±1.0 GeV ll26l . The fitted spectra are shown in Figure |9] We estimate the 
systematic uncertainty by vaiying M o within this range and repeating the fit, it amounts to 

A) 

±2 GeV on the determination of M v ± . 

Xi 

The relative mass accuracy for the two models scales with production cross section as ex- 
pected. We notice that the deterioration of the mass accuracy due to BS is larger for model I 
(+50.4%), where the X\ mass is significantly smaller compared to the beam energy, compared 
to model II (+28.4%), since the X\X\ threshold sits close to y/s and the BS effects are dumped 
by the fall of the production cross section near threshold. Then, we repeat the fit accounting 
for both the direct chargino production contribution to the W energy spectrum and that from 
— > % ± v e decays. These decays still offer some sensitivity to the chargino mass but cannot be 
distinguished from chargino pair production in the event selection. At this specific benchmark 
point the selectron left mass M eL cannot be measured with the electron energy spectrum in the 
e + Xi e Xi final state since the spectrum is dominated by the SrSr process. It could be measured 
using the mixed e + X[h°X[h° and possibly e^fW^V^j modes at 3 TeV and in a threshold scan. 
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Table 5: Statistical accuracy on %\ mass from a fit to the E w ± endpoints for 2 W + missing 
energy inclusive SUSY events and 2 ab _1 of integrated luminosity at 3 TeV, under 
different assumptions 



Particle 


Mass 


No Rad 


ISR 


ISR+BS 


ISR+BS 


ISR+BS 




(GeV) 






SE/E=0 


=0.025 


=0.050 


Model I 






2-par Fit 














xt 


643.2 


± 1.4 


± 1.7 


±2.6 


±3.8 


±4.1 


3-par Fit 














xt 


643.2 


± 2.0 


± 2.3 


±3.5 


±4.8 


±5.2 




1102.2 


±48 


±50 


±56 


±63 


±66 



We perform fits where Mg L is either kept at its model value or it is treated as a free parameter, 
while the xt mass and the relative contribution of the two processes to the W spectrum are kept 
free. Results are summarised in Tableland Figure ITOl 

3.3.2 X2 -> h°Xi 

The x\ ~^ h X\ decay of the chargino pair yields two Higgs bosons plus missing energy. The 
same final states originates also from e + e~ — > ViVi — > X^ytX^yt- F° r the parameters of model I 
My e = 1097.2 GeV, the e + e~ — > V/Vg production cross section is 14.5 fb. The sneutrino channel 

Table 6: Statistical accuracy on x% mass from a fit to the E h o endpoints for X2X1 signal events 
and 2 ab _1 of integrated luminosity at 3 TeV, under different assumptions 



Particle 


Mass 
(GeV) 


No Rad ISR ISR+BS ISR+BS ISR+BS 
8E/E=0 =0.025 =0.05 


Model I 

X° 2 


643.2 


±1.01 ±1.17 ±2.58 ±3.59 ±4.54 


Model II 

X° 2 


1064.2 


±10.64 ±11.12 ±16.71 ±19.04 ±23.42 



accounts for 68% of the 4-jet, hh + missing energy inclusive SUSY sample. In model II, the 
tree-level sneutrino production cross section is 6 fb which drops to 2.3 fb accounting for ISR 
and BS and the signal purity is 89% . This final state is expected to be virtually immune from 
irreducible SM backgrounds, since the cross section for double WW fusion light Higgs produc- 
tion, e + e~ ->■ h°h a v e v e is only ~1 fb for 1 15 < M h < 130 GeV [27 ]. Fit results on signal events 
are summarised in Tableland the fitted spectra given in Figure [TT] Again, the relative increase 
of the statistical uncertainty due to BS is larger for model I compared to model II. We repeat the 
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°0~ U 200 ' 400 ' 600 ' 800 ' 1000 ' I2OO °0^ 200 ' 400 ' 600 ' 800 ' 1000 ' ^200 

E h (GeV) E h (GeV) 

(a) ISR Only (b) ISR + BS 

Figure 1 1 : h° energy spectrum in the 4-jet, 2 h + missing energy final state for X2X2 S1 g na l events 
in model I at 3 TeV with only ISR (left) and ISR+BS (right). The points with error 
bars are the simulation and the continuous line the fitted spectrum. 

Table 7: Statistical accuracy on %® mass from a fit to the E w ± endpoints for 2 h + missing energy 
inclusive SUSY events and 2 ab _1 of integrated luminosity at 3 TeV, under different 
assumptions 



Particle 


Mass 
(GeV) 


NoRad ISR ISR+BS ISR+BS ISR+BS 
8E/E=0 =0.025 =0.050 


Model I 

X°2 

ve 


643.2 
1097.2 


± 2.5 ± 3.2 ± 5.3 ± 6.8 ± 8.4 
±43.4 ±51.8 ±52.7 ±60.4 ±69.3 



fit accounting for both the direct x\ production and the irreducible Va SUSY background con- 
tribution to the W energy spectrum for model I where its contribution is more important. Since 
the sneutrino decays through a the h° energy spectrum retains some sensitivity to the X2 
mass. A multi-parameter fit with the %® and v<? masses treated as free parameters, together with 
the relative contribution of the two processes to the h° spectrum is performed (see Figure [12}. 
Results are summarised in Table [7] 

3.3.3 xi -> W ± x°, xt -)• W ± x!, %t -> fPZi 

The X2 chargino has one of largest production cross section for supersymmetric particles. Con- 
trary to the lighter states discussed above, there is no dominant decay channel and the analysis 
has to include several final states. We consider here W ± , Z° and h° spectra from combina- 
tions of decays with with both 6- and 8-parton final states. The channels we have consid- 
ered are not exhaustive but are representatives of the topologies offered by decays of X2X2 
pairs. The 6-jet topology accounts for 31% of the total yield of X2X2 P arrs - The W + W~h° 
+ missing energy final state receives two contributions from signal X2X2 P a * r production: 
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°0 -J 200 ' 400 ' 600 ' 800 ' 1000 ' ^1200 

E„ (GeV) 



pi i B i 
200 
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E„ (GeV) 



(a) ISR Only (b) ISR + BS 

Figure 12: h° energy spectrum in the 4-jet, 2 h + missing energy final state for inclusive SUSY 
events in model I at 3 TeV with only ISR (left) and ISR+BS (right). The points with 
error bars are the simulation, the continuous line the fitted spectrum and the dashed 
line the fitted contribution from X1X2 events. 

Table 8: Statistical accuracy on xf mass from the combination of fits to the E w ± and E h o spec- 
tra for X2X2 S1 g na l events in 6-jet WWh and 8-jet WWZh topologies and 2 ab _1 of 
integrated luminosity at 3 TeV, under different assumptions 



Particle 


Mass 


No Rad 


ISR 


ISR+BS 


ISR+BS 


ISR+BS 




(GeV) 






8E/E=0 


=0.025 


=0.050 


Model I 














xf 


916.7 


±2.2 


±2.8 


±3.6 


±3.9 


±4.2 



e + e~ -)• X2X2 ~> W + X% W~Xi ~> W + h°XiW~Xi with h -)■ bb, W -)• qq' and e + e~ -)■ X2X2 ~> 
h°XiW~Xi h°W + XiW~X\ with h — > bb, W — > qq'. Each event has two W and one Higgs 
boson. The purity in X2X2 P an ' s i s 77 % w hh other SUSY contribution to this final state com- 
ing from X3X4 an d X\Xa ■ The W spectrum combines Ws produced in the direct x 2 decay 
into two different states, X2 an ^ X\ • m the final state also the Higgs boson is sensitive to the 
X 2 mass, through the contribution from the x 2 ~^ h°X\ decays, while part of the bosons come 
from the x\ decay. In both cases we perform a 2-par fit leaving both the M^± and the fraction 
of the two contributions free (see Figure [13). Then, we consider processes with four bosons 
yielding the 8-jet topology, which accounts for 19% of the X2X2 h na l sta tes. The Z°h°W + W~ 
final state has 85% contribution from X2X2 P arrs with the remaining mostly due to Xi%a- We 
use both the Z° and the h spectrum to perform the mass fits. In this topology the Z energy 
spectrum from the process e + e~ — > X2X2 ~~ ^ X^X\^~ X2 ~~ ^ Z°W + W~h°X[Xi receives almost 
exclusively contributions from the direct x 2 ~^ decay. For the fit to the h° spectrum we 
include the contribution from the xf ~^ h°xt an d that from x 2 ~^ h°Xi> leaving the contribution 
of the two processes free, as done for the 6-jet topology (see Figure [13]). Results are combined 
and the statistical uncertainties are summarised in Table [8] 
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3.3.4 xg^W^tf 

The e + e~ — > X3X4 process has a sizable cross section in model I and a good number of signal 
events can be reconstructed. The decay channel considered here is %® 4 — > W ± xT — > W ± W^fXi 
which gives an 8-jet topology with, 4 W bosons and missing energy. This final state is chal- 
lenging, due to its large jet multiplicity, but the signature is striking and there are essentially no 
SUSY or SM irreducible background processes contributing to it. In each event two of the W are 
sensitive to the X34 mass and the others to the X\ mass. Since this can be precisely measured 
in the 2 W + missing energy channel, it is safe to assume here that its mass is known. Results 
are given in Table [9] and Figure [14] The structure of the cMSSM implies that the two heaviest 
neutralinos, %3 and X4 w& nearly degenerate. However, this is not a general feature of super- 
symmetry and it does not apply to less constrained scenarios, such as the MSSM, where the %® 
- %® mass splitting can be ~20% of their masses, or larger. Therefore, we repeat the fit, leaving 
the masses of X3 an d X4 f ree an d test the result for a mass splitting of 10 GeV, as in our model I, 
and one of 40 GeV obtained by decreasing the #3 mass. We find a resolution on the X3 - X4 
mass splitting of ~ 15-18 GeV and masses which are apart ~35 GeV and more can be resolved. 




E w (GeV) E„ (GeV) 

(a) 6-jet WWh (b) 6-jet WWh 




E z (GeV) E„ (GeV) 

(c) 8-jet WWZh (d) 8-jet WWZh 

Figure 13: Energy spectra for X2X2 S1 g na l events: (upper row) in the six jet WWh final state and 
(lower row) in the eight jet WWZh final state with ISR and BS effects for model I at 
3 TeV The points with error bars are the simulation and the line the fitted spectrum. 
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(a) ISR Only 



1000 1200 
E w (GeV) 




1000 1200 
E„ (GeV) 



(b) ISR + BS 



Figure 14: energy spectrum in the 8-jet topology, 4W + missing energy final state at 3 TeV 
with only ISR (left) and ISR+BS (right) for X3X4 signal events in model I. The points 
with error bars are the simulation and the line the fitted spectrum. 



Table 9: Statistical accuracy on X34 mass from a fit to the E w ± endpoints for X3X4 signal events 
and 2 ab -1 of integrated luminosity at 3 TeV, under different assumptions 



Particle 


Mass 
(GeV) 


NoRad ISR ISR+BS ISR+BS 
SE/E=0 =0.05 


Model I 
v o 

A,4 


905.5 
916.7 


±7.1 ±7.9 ±12.6 ±15.1 
±8.7 ±8.9 ±13.0 ±15.4 



3.4 Energy Resolution Effects 

The effect of the jet energy resolution on the X\ masses for both models are shown in Figure [T5l 
in terms of the relative change of the statistical accuracy of the mass determination as a function 
of the parton energy resolution 8E/E. As expected, the effect is larger when beam radiation 
is not considered. We establish a quantitative criterion for the energy resolution requiring that 
the contribution from the parton energy resolution to the statistical mass accuracy should not 
exceed the combined effect of ISR and beamstrahlung. We assume that these contributions adds 
quadratically and derive the limits to 8E/E for the different channels which are summarised 
in Figure [16] Then, we consider the effect of the parton energy resolution on the di-jet 
invariant mass. We study the specific case of the Z° identification in the 8-jet e + e~ — > X2X2 ~* 
Z°x^W~ Xi — * Z°W + X\W~ X\'> Z ~^ <7<7> W — > qq' process, which is the process most sensitive 
to di-jet mass resolution in model I, due to the large W yield. We identify the Z° boson requiring 
that the di-jet mass is compatible with the nominal mass, Mz=91.2 GeV, within 3 Oz, where 
Gz is the peak Gaussian width measured on genuine Z° — > qq di-jets. We vary the Gaussian 
parton energy resolution (see Figure [T7J and study the purity in real Z° bosons selected by this 
selection. Since W ± bosons are dominant in 8-jet topology SUSY events (53%) and Z° bosons 
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model 
• No Radiation 
■ ISR only 
- ISR + BS 




model II 
• No Radiation 
■ ISR only 
' ISR + BS 




.02 0.04 0.06 0.08 0.1 0.12 



(a) Model I 



(b) Model II 



Figure 15: Relative change of the statistical accuracy on the %t mass as a function of the jet 
energy resolution from one-parameter % 2 fits to the energy spectrum for signal events 
in model I (a) and II (b). 



0.15 

SE/E 



Figure 16: Values of the relative energy resolution 8E/E at which its contribution on the mass 
statistical uncertainties equals that of ISR and beamstrahlung for the channels con- 
sidered in this study. Model I is shown by the open histogram and model II by the 
grey histogram. 




(a) <5£/£=0.025 





(b) S£/£=0.050 



(c) SE/E=0.015 



Figure 17: Di-jet invariant mass in 8-jet inclusive SUSY events in model I for various values of 
parton energy resolution. 
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Table 10: Gaussian width of the Z peak and purity of Z decays selected from their compatibil- 
ity with the Z° mass in 8-jet inclusive SUSY events in model I for various values of 
parton energy resolution. 



8E/E 


a z 


Purity 




(GeV) 




0. 


1.64 


0.877 


0.025 


2.42 


0.720 


0.040 


3.33 


0.418 


0.050 


AAA 


0.290 


0.075 


5.34 


0.220 



make only 15.5% of the di-jets, the leakage from the W ± peak is important already for moderate 
values of the energy resolution, as shown in Table |T0j 

Operating the collider with polarised beams may be important to improve the statistical accu- 
racy in the determination of heavier states such as %® and %®. However, given the broad scope of 
the research program at a multi-TeV collider various states of polarisation will be likely selected, 
making our assumption a fair estimate for the total cumulative statistics of signal events. 

3.4.1 Validation with Full Simulation and Reconstruction 

In order to validate the results obtained above at generator level, accounting only for a simple 
energy smearing, the xt'X.i anc ^ %2%2 analyses are repeated on fully simulated and reconstructed 
events to verify the accuracy when accounting for the reconstruction effects in full. 

A sample of inclusive SUSY events for model I, generated with ISR+BS, corresponding to 
0.5 ab~' of integrated luminosity, is fully simulated and reconstructed using the CLIC version of 
the ILD detector concept. First, the 4-jet + missing energy events are reconstructed. Events are 
pre-selected requiring a visible energy 250 < E tot < 1 800 GeV, an energy in charged particles 
larger than 150 GeV, transverse energy larger than 200 GeV, a jet multiplicity 2< Nj ets <5 and at 
least 20 charged reconstructed particles. Jets clustering is performed using the Durham jet algo- 
rithm ESI , with y cut = 0.0025, on the reconstructed particle flow objects of the Pandora particle 
flow package ll29l . These events are then forced into four jets and the di-jet invariant mass for all 
the three possible pairings is computed. The jet pairing minimising the difference between the 
di-jet invariant masses is selected, provided the mass difference is below 20 GeV. The resulting 
mass distribution on generated SUSY di-boson events is shown in Figure [18] The fraction of 
W + W~ , Z°Z° and h°h° events is extracted by a % 2 fit to the di-jet mass distribution. The W and 
Z° mass peaks are parametrised as Breit-Wigner functions convoluted with a Gaussian term de- 
scribing the experimental resolution. The mass and width values of the Breit-Wigner functions 
are fixed to their generated values, while the total area and the width of the Gaussian resolution 
terms are left free in the fit. The h peak, which has negligible natural width, is modelled as the 
sum of two Gaussian curves, one representing the correctly reconstructed signal events, centred 
at the nominal Mf, value, the second describing decays where the mass has a lower reconstructed 
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" " "40 60 80 100 120 140 160 180 200 

Figure 18: Di-jet invariant mass for di-boson inclusive SUSY fully simulated and reconstructed 
events for model I with the fitted contributions from W , Z° and h 

Table 11: Fraction of W ± , Z° and h° bosons from the fit to the di-jet invariant mass distribution 
in 4-jet + missing energy inclusive SUSY events for model I for 0.5 ab _1 of integrated 
luminosity, compared to the generated values. 



Boson 


Fitted 


Simulated 




Fraction of Evts. 


Fraction of Evts 


W ± 


0.650 ±0.011 


0.645 ± 0.005 


z° 


0.040 ± 0.009 


0.020 ± 0.002 


h° 


0.215 ± 0.010 


0.243 ± 0.003 



value due to semi-leptonic b decays. The central value, width and fraction of events in this sec- 
ond Gaussian is extracted by a fit to a pure sample of decays into h bosons and fixed in the fit, 
while the Gaussian width of the main peak is kept free. Results are given in Table [TT] and the 
fitted functions are overlayed to the reconstructed spectrum in Figure [T8l 

Then, events with di-jets compatible with the WW hypothesis are selected. The total selec- 
tion efficiency is 60% for X\X\ ~^ W + X\W~ X\ —> <7<7 'Xiffl 'X\ signal events. This efficiency 
is independent on the W energy and the reconstruction and selection criteria do not introduce 
any significant bias to this distribution (see Figure H9l . The sample of selected 4-jet WW can- 
didate events has a purity of 86% and consists of 78% signal X\X\ ~^ W + X\W~Xi> 7.5% 
e^e^ ->■ X\X\ v e% and 0.2% of other 4-jet WW SUSY processes (see the left panel of Fig- 
ure l20l . The fraction of background fake-WW events is obtained from the fit to the di-jet in- 
variant mass distribution and the shape of their di-jet energy spectrum directly extracted from 
the reconstructed events, using the di-jet mass side-bands 40 < Ejj < 60 GeV and 140 < Ejj < 
160 GeV, and subtracted. 

The multi-parameter fit to the X\ an d masses is repeated on the background-subtracted 
W^ 1 energy distribution of selected fully simulated and reconstructed events (see Figure|2T}. The 
result is M v ± = (643 ± 14) GeV and M ; ± = (1 100 ± 104) GeV, where the statistical accuracies 

A1 " I 
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Figure 19: Efficiency for event reconstruction and selection on fully simulated and reconstructed 
4-jet Xi%i ~^ W + X®W ~X\ — > qq'XiW'Xi signal events as a function of the W en- 
ergy. 



600 




Mjj (GeV) Mjj (GeV) 

Figure 20: Di-jet invariant mass for selected inclusive SUSY fully simulated and reconstructed 
events for model I (left) with the different sources highlighted and (right) with the 
background and signal components fitted. 

are consistent to those obtained on a smeared simulated spectrum of equal statistics, M ± = (643 
± 12) GeV and M s ± = (1 100 ± 1 10) GeV. Finally, we select events with di-jets compatible with 
the hh hypothesis. The total selection efficiency is 69% for X2X2 ~~ ^ h°Xih°Xi -> bbx\bbx\ The 
multi-parameter fit to the X2 an d Vl masses is repeated on the background-subtracted h energy 
distribution of selected fully simulated and reconstructed events (see Figure 1221 . The result 
is M v o = (643 ± 26) GeV and My, = (1097 ± 148) GeV, where the statistical accuracies are 
consistent to those obtained on a smeared simulated spectrum of equal statistics, M o = (643 ± 
21) GeV and M Vf = (1097 ± 123) GeV. The degradation of the statistical accuracy in the analysis 
of the fully simulated and reconstructed data compared to that on the smeared generator-level 
events, is likely due to the larger background from fake 4-jet hh events, which is not included 
at generator level. In this study we do not consider jet flavour tagging for consistency with the 
simple procedure adopt for the generator level study. However, by applying ^-tagging to the four 
jets, the background, which does not contain two light Higgs bosons, can be largely reduced. 
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Figure 21: Fits of fully simulated and reconstructed events in model I: (left) W ± energy spec- 
trum of selected 4-jet WW candidate events energy spectrum and (right) background 
subtracted spectrum (points with error bars) with the result of the 3-par fit (continuous 
line). 



rzh v, 

| | other 4-j hh SUSY 




200 400 600 800 1000 1200 

E h (GeV) 




Figure 22: Fits of fully simulated and reconstructed events in model I: (left) h° energy spectrum 
of selected 4-jet hh candidate events energy spectrum and (right) background sub- 
tracted spectrum (points with error bars) with the result of the 3-par fit (continuous 
line). 



In conclusion, the analysis of fully simulated and reconstructed SUSY events where 4-jet, 
WW and hh candidates are selected based on the reconstructed topology and di-jet invariant 
mass shows that reconstruction efficiencies are quite large and flat with the boson energy. Back- 
grounds from other final states can be reliably estimated and subtracted in a model-independent 
way. The statistical accuracies obtained on the extraction of gaugino and slepton masses from 
these data are found to be comparable to those from the smeared generator, once the reconstruc- 
tion efficiencies are taken into account. 
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Figure 23: Purity for the X\X\ an( ^ X2X2 P rocess in 4-jet final state as a function of the y/s 
energy. 

4 Mass Determination by Threshold Energy Scans 

An e + e~ linear collider with tunable beam energy can determine the sparticle masses by per- 
forming energy scans of their pair production cross section near threshold. In principle, this 
method provides a better mass accuracy, compared to the kinematic end-point method discussed 
above. Threshold energy scans put significant requirements on the machine performance and 
versatility. Not only the beam energy needs to be varied over a broad range, but, since the cross 
section at threshold is small, a large luminosity must be preserved in lower energy operation. 
Beamstrahlung effects are important at threshold, while SUSY background are reduced, at least 
for the lighter states. We study the processes 

• e+e- -> xtXi -> W+xy-X®; W -> qq', 

• e + e~ — > X2X2 ~+ h°Xih°Xi', h — > bb, 

• e + e- -> xtxt -> W+x^W'Xi -»• W + h°ZiW~Xi> h^bb,W -> qq' 

• e+e- -> £ 2 + * 2 + -> fc^+W-^ -> ^W+^fW-^ ; h -)■ &5, W -> <?<f 
for model I. For model II we study 

• e+e- -> -»• W^W"*?; W -> <?<?' 

We assume a total integrated luminosity of 3 ab _1 , where the 2 ab -1 taken at the maximum 
energy, as assumed above, are supplemented by 1 ab -1 of statistics dedicated to the scan of 
sparticle pair production thresholds at lower energies. At, or below, y/s = 2 TeV, the X\X\ 
and #2 #2 P a * r production saturates the final states with 4-jet and WW or hh, respectively, since 
slepton production is below threshold due to the larger &l and Vg masses, as shown in Figure 1231 
A 1-par x 2 fit to the cross section values at the chosen operating energies is performed to extract 
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Fitted Mass (GeV) Fit Stat. Error (GeV) 



Figure 24: Toy test validation of the threshold scan fit results for %f for model I. The corre- 
sponding x 2 At result is (643.2± 0.68) GeV. 
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Figure 25: Sensitivity -^j^m to sparticle masses in the threshold scan as a function of yfs for 
model I. 



the mass value. The mass and its uncertainty are obtained by assuming a given number of 
cross section measurements of the relevant pair production process at the y/s values for the scan 
points. The cross section is computed both at Born level and also adding ISR and beamstrahlung 
effects using Pythia. The fit results are validated using toy tests. In these we repeat the fits 
by varying the cross section within its statistical uncertainty at each ^/s value and we plot the 
result. We verify that the result is centred on the simulated mass and its width is consistent 
with the x 2 fit uncertainty (see Figure [24). In order to define a suitable scan strategy, we first 
study the sensitivity to the sparticle masses as a function of ^fs. The sensitivity is denned as 
1 / y/odo I dm, where a is the pair production cross section and m the mass, as in ref. QUI . We 
compute the cross section a at various -y/s values for a set of closely spaced masses and obtain 
the derivative do /dm of the change of the cross section at each energy per unit of mass change. 
Results are shown in Figure |25j which indicate that the maximum of the sensitivity to the mass 
is achieved near threshold. The number of scan points and the share of the statistics among them 
is optimised by studying the mass uncertainty obtained from the fit for different assumptions. 
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Table 12: Statistical accuracy on sparticle masses from energy scans under different assumptions 
for 2 ab _1 of integrated luminosity at 3 TeV and 1 ab _1 at two energies near threshold. 



Particle 


Mass 
(GeV) 


Born 


ISR 


ISR+BS 


ISR+BS 
+Bkg 


w/Pol 
(+0.8/0) 


w/Pol 
(+0.8/-0.6) 


Model I 
















xt 


643.2 


±0.6 


±0.6 


±0.7 


±0.7 


±0.5 


±0.4 


X°2 


643.1 


±4.3 


±13.8 


±24.1 


±25.6 


±23.9 


±18.1 


xt 


916.7 


±0.8 


±0.9 


± 1.3 


± 1.4 


± 1.1 


±0.9 


Model II 

xt 


1062.2 


± 6.2 


±6.4 


±6.9 




±5.1 


±2.8 



We find that it is preferable to concentrate the luminosity in a small number of scan points. For 
example, the statistical accuracy on the mass of the %f~ in the model I varies from ±0.85 GeV, 
obtained for a four-point scan (13 10< y/s <1950 GeV), to ±0.45 GeV, when the luminosity 
is split between just two points, one of which at the peak of the sensitivity (y / ?=1350 GeV) 
and the second close to threshold (^=1310 GeV). This confirms the findings of [30] for lower 
sparticle masses and different luminosity spectrum. Finally, we consider the option of operating 
the collider with polarised beams. Results are summarised in Table [T2l In all cases, except the 
the mass accuracies obtained with a dedicated threshold scan improve on those resulting 
from the kinematic edge analysis at 3 TeV by factors of 2 or more. The use of polarised beam 
further improves these accuracies, effectively compensating for the loss of sensitivity due to ISR 
and BS. 

5 Conclusions 

The determination of chargino and neutralino masses in high-mass SUSY scenarios with two- 
body decays into W , Z° and h bosons provides us with a mean to quantify the effect of ra- 
diation, by ISR and beamstrahlung, and parton energy resolution on the accuracy achievable in 
supersymmetric particle mass measurements at a multi-TeV e + e~ linear collider. In our anal- 
ysis both fits to boson energy spectra and threshold scans are considered for fully hadronic 
final states. Results from generator-level quantities are validated using fully simulated and re- 
constructed events in the IV + ^missing and h°h Q + ^missing final states. Not accounting for 
reconstruction efficiencies, estimated to be ~60% in four jet final states, the mass of charginos 
and neutralinos can be determined from the kinematic edges of the boson energy in inclusive 
SUSY event samples to a relative accuracy in the range 0.3% to 1.0% (0.6% - 1.0%) in absence 
of radiation and energy resolution effects to 0.8% to 1.7% (1.1% - 2.0%) accounting for ISR, 
BS and realistic energy resolution for the benchmark with particle masses in the range 600 - 
900 GeV (>1000 GeV), respectively, with 2 ab _1 of integrated luminosity at y/s = 3 TeV. The 
relative increase of the statistical uncertainty of the mass measurement is larger for the model I 
which has the sparticles masses far way from pair the production thresholds. However, in abso- 
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lute terms the larger production cross sections in this model yield better statistical accuracy in 
the mass determination. By adopting the criterion that the degradation to the mass measurement 
statistical accuracy from the partem energy resolution should not exceed that induced by ISR 
and BS, we derive the requirement of a relative energy resolution for jets, 8E/E <0.05. If the 
accelerator can operate at energies below the nominal yfs (down to ^=1310 GeV for model I 
and -^=2200 GeV for model II) with comparable performance to collect about one third of the 
statistics at centre-of-mass energies close to the kinematic thresholds for sparticle pair produc- 
tion, the mass accuracies from these threshold scans improves by factors of 2 or more compared 
to those obtained from study of the kinematic edges at the maximum y/s energy The availability 
of polarised beam in the scan further improves these accuracies, effectively compensating for the 
loss of sensitivity due to the effect of ISR and beamstrahlung. 
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